CERN-TH/97-361 
hep-ph/9712343 
November 1997 

The Role of Scalar and Pseudoscalar Fields in Determining 
Nucleosynthesis Bounds on the Scale of Supersymmetry Breaking 



Tony Gherghetta* 
00 ' Theory Division, CERN, CH-1211 Geneve 23, Switzerland 

o\' 

^ ■ Abstract 

The effect of spin-0 goldstino superpartners is considered on the nucleosynthe- 
sis bounds arising when a superUght gravitino appears as an effective mass- 
^N \ less neutrino species. When the scalar and pseudoscalar superpartners are 

^f^ \ relativistic they will decouple at much later times than the goldstino and 

''^ ■ consequently will be the dominant effect when obtaining a nucleosynthesis 

bound on the scale of supersymmetry breaking. Assuming that the scalar 
and pseudoscalar fields decouple at a temperature no later than O(IOO) MeV, 
L-. \ then typically the scale of supersymmetry breaking yF ^ 60 TeV. This cor- 

ix. 



en 



X 



responds to a lower bound on the gravitino mass w.3/2 ^ 1 eV. 
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I. INTRODUCTION 

If supersymmetry is realised in nature then it is presumably broken at some energy scale 
and subsequently transmitted indirectly to the low energy observable mass spectrum [|l|. 
While naturalness arguments essentially constrain the amount of supersymmetry breaking 
in the low energy mass spectrum, the scale at which supersymmetry is broken remains largely 
undetermined. In general, realistic supersymmetry breaking scenarios can be parametrised 
by a spurion chiral supermultiplet (0, ip, F) in which a nonzero vacuum expectation value for 
the F-term gives rise to the spontaneous breaking of supersymmetry at the scale yF. The 
fermion, ip represents the massless Nambu-Goldstone fermion which becomes part of the 
massive gravitino after the superHiggs effect. This fermion is referred to as the goldstino G. 
The gravitino mass 777,3/2 is determined purely by gravitational interactions and is fixed by the 
relation F = A/37773/2MP, where the reduced Planck mass Mp = a/SttGtv = 2.4 x lO^^GeV 
is related to the gravitational Newton constant G^. When the gravitino mass 7713/2 — *■ 0, 
the goldstino becomes the dominant component of the gravitino during interactions. 

The mass of the spin-0 superpartner 0, is however not fixed and can in general depend on 
nonminimal terms in the Kahler potential as well as terms in the superpotential. When the 
complex scalar is decomposed into real fields (p = {S+iP)/\/2, we need to consider the scalar 
mass ms and pseudoscalar mass Trip as additional parameters to the supersymmetry breaking 
scale F. In this framework we will be interested in placing a lower bound on the scale of 
supersymmetry breaking from the effect of a superlight gravitino together with its scalar 
superpartners contributing as extra effective neutrino species during the nucleosynthesis 
epoch. 

The effect of a superlight gravitino in the early universe was considered in a previous 
calculation |0. There it was assumed that the spin-0 superpartners were essentially mass- 
less [ms, mp -C T, where T is the temperature during the nucleosynthesis epoch) and that 
they decouple before the superlight gravitino. More recently, there have appeared additional 
papers [^]-|^ in which the general question of determining an effective gravitino Lagrangian 
have been addressed. In addition these references investigated the case of massive super- 
partner masses ms, m,p ^ T and noticed qualitatively different behaviour for scattering 
cross sections than the superlight scenario. However they did not address the question of 
superlight scalar superpartners during nucleosynthesis. 

In this work we will consider the effect of the scalar superpartners in determining the 
nucleosynthesis bounds on the supersymmetry breaking scale F and in particular consider 
the case of superlight masses for the scalars. It will be shown that S and P decouple 
much later than the goldstinos, so that the scalar and pseudoscalar contributions actually 
strengthen the bound considered in [Q], when ms and irtp are superlight as well. While 
we provide a more complete study of the effect of the scalar and pseudoscalar fields during 
nucleosynthesis, the result for superlight S,P masses will agree with 0. In section 2 we 
will consider the scattering cross sections for S, P and G which are relevant during the 
nucleosynthesis epoch. This will then allow us, in section 3, to obtain an expression for 
the corresponding decoupling temperature. We will see that this leads to a bound on the 
supersymmetry breaking scale, F and consequently a lower bound on the gravitino mass 
7773/2. The conclusion and final comments appear in section 4. 



II. SCATTERING CROSS SECTIONS 

We will assume that all gravitino interactions are dominated by a superlight goldstino in 
determining the gravitino scattering cross sections. The associated scattering cross sections 
of the scalar superpartners will also be considered when they are superlight as well. We 
will make the usual assumption that a massless and weakly interacting particle decouples 
when the relic abundance is frozen. Thus it will be sufficient to consider processes that 
change particle number. The scattering cross sections can be obtained once all the effective 
interactions between the relevant particle species is known. While the full supergravity 
formalism was employed in 0,0], it will be simpler instead to use the effective Lagrangian 
approach adopted by [^ where the effects of gravity are formally decoupled in the limit 
Mp — > oo while keeping F fixed Q. During the nucleosynthesis era the particle species 
that are assumed to be in thermal equilibrium besides 5, P and G will be the photon (7), 
three families of massless neutrinos (z/j) and fermions, including at least electrons (e^) and 
possibly muons (yU^). The temperature of the universe at these times will be in the range 
Oil - 100) MeV. 

A. Goldstino annihilation 

The cross section for goldstino annihilation in the massless S*, P limit was previously 
calculated for various channels in 0J^. For goldstino annihilation into photons the result is 

1 M^ 
aiGG -^ 77) = -^^s\ (1) 

^ "^ 647r F4 ' ^ ^ 

where s denotes the Mandelstam variable for the total energy in the center of mass frame and 
M generically denotes the gaugino mass (the effects of neutralino mixing will be ignored in 
this paper). The result (^ is valid in the limit 711^/2, rns,rnp -C y/s <^ M. It is also 
straightforward to obtain the annihilation cross section in the limit ^/s <^ ins.,mp,M for 
which one finds 

aiGG -^ 77) = -^^. (2) 

This expression for the annihilation cross section agrees with the result derived by [0,0, 
where it was noted that the cancellations of the leading order terms change the massless 
S", P limit cross section (|I]) . These cancellations are obviously absent when one derives the 
cross section in the massless S", P limit and there is no inconsistency with the previous 
literature [0,|],01, in which the massless limit was always the underlying assumption. 

Similarly, the Goldstino annihilation into fermions was calculated in p|,^. As shown in 
^ the cross section for annihilation into fermions is ambiguous up to an unknown parameter 



^An alternative approach using a nonlinear realisation of the supersymmetry algebra can also be 
used U. 



which represents the neglect of higher derivative terms in the effective Lagrangian approach. 
While this ambiguity will be seen to have no direct consequences on the results obtained in 
this paper, we will nevertheless use a minimum value of the cross section, namely H] 
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in order to derive a model independent bound. 

The above annihilation channels into photons and fermions are not the only possible 
channels. Since the goldstino superpartners 5* and P are also assumed to be in thermal 
equilibrium we need to also consider the annihilation channels GG -^ SS, PP, SP. The 
effective Langrangian governing these interactions is given by 
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This leads to the following Goldstino annihilation cross sections into the scalar superpartners 
S, P in the limit m^, mp -C ^i 
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Clearly in the massless limit these cross sections are significantly smaller than ([l|) and 
will contribute negligibly to the overall goldstino annihilation rate. 



B. S and P annihilation 



The annihilation channels available for S and P are similar to those for the goldstino, 
except that unlike goldstinos there are no fermion couplings in the leading order 1/ F . Thus 
the effective Lagrangian describing their interactions will be the goldstino Lagrangian (^ 
together with the interaction with photons H 
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In the limit that ms,Tnp <^ ^/s this leads to an annihilation cross section 

5 IMI^ 
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obtained via t and u-channel photon exchange. This cross section agrees with the result 
obtained by employing the full supergravity formalism [§]. Notice that the dependence on 
the Mandelstam variable s is qualitatively different from that of (|1|) and we can already 
see that 5* and P interact more strongly with the thermal plasma when T < M. This will 
ultimately correspond to S and P decoupling later than the goldstino. 

There are also possible annihilation channels into goldstino pairs due to the interactions 
(^). These are the inverse of the processes considered in (^^ and give rise to similar 
cross sections for annihilation into goldstino pairs, namely a{SS -^ GG) = Aa{GG — »• 
SS) and similarly for P. The overall factor of 4 comes from the averaging over initial 
polarisations of the goldstino in the expressions for the unpolarised cross section (^^. 
Again this annihilation channel will have a negligible effect on determining the decoupling 
temperature. Other channels such as SG — > PG will also not be important. 

The form of the S, P couplings to photons (|-|^) also allows for the channel Sf^jf. A 
photon is exchanged in the t-channel and couples to fermion pairs in the standard way. In 
the limit ^/s ^ nis, mj where rrif is the fermion mass one finds that 
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where a is the electromagnetic fine structure constant and / is a Dirac fermion (similarly 
for P with the substitution ms —>■ mp). Notice that the dependence on the center of mass 
energy ^/s is now logarithmic and no longer a power law like previous annihilation channels. 
This channel (known as the Primakoff reaction) was also considered in [^] but in the limit 
rrif ^ E^ ^ Es and is also important for establishing bounds on the gravitino mass from 
supernova cooling |10,B|. When ms — > the cross section becomes singular due to collinear 
singularities. In order to regulate this infrared divergence we can introduce a thermal mass 
M^ for the photon in the background plasma. Without going into a full finite temperature 



calculation this will be a good approximation to leading log order |]ri| . Thus, in the massless 



S limit the cross section (|TT|) becomes 
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and similarly for P. This cross section will become important in obtaining the decoupling 
temperature of the S", P superpartners. 

Similarly S and P can both scatter off the background photons via a t and s-channel 
photon exchange with various particles in the final state. The scattering cross section for a 
photon in the final state in the massless 5, P limit is 



cr(S'7 — » 7P) = cr(P7 — > 76*) 
Similarly the cross section for pair production is 
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The process (|T3D is comparable to (0), while ([M]) is subdominant compared to (jT^) and 
neither will have a large effect in the determination of the decoupling temperature. 

Up to now we have implicitly neglected the possibility of cubic scalar couplings in the 
previous goldstino and scalar channels. In general these cubic scalar interactions arise from 
nonminimal terms in the Kahler potential and superpotential terms ||^. While in all realistic 
models the cubic scalar couplings vanish it is interesting to comment on their effect if the 
scalar annihilation channel were to include such a term. For example, the effect of the scalar 
cubic couplings on the cross section (|T0|) can be seen if we parametrise the interaction as 
C = —As/3\ S^. Including this term allows for an s-channel S exchange diagram contributing 
to the annihilation into photons. In the massless S limit we obtain 

The interesting feature of the above cross section is that while at high energies we repro- 
duce the result (ll^), at low energies the cross section may grow due to the 1/s behaviour. 
Cosmologically this would mean that the scalar S would actually come back into thermal 
equilibrium at late times if the coupling As were large enough. Similarly a cubic scalar 
coupling of the form ~ S^P would modify (0) to a form like (|15D via a u-channel scalar 
exchange diagram. 

There are also dimension-6 terms ~ SSF^^F^" (similarly SP, PP) which in a more gen- 
eral framework must also be included^ The coupling is controlled by the second derivative 
of the gauge kinetic function. Again for simplicity we will assume this term to be absent. 



III. THE DECOUPLING TEMPERATURE 

The decoupling temperature of a particle species in the early universe is obtained in 
general by solving a set of coupled Boltzmann equations. However, in practice a reasonable 
estimate of the decoupling temperature can be obtained by checking when the interaction 
rate of a particular particle species begins to fall behind the Hubble expansion rate H . The 
interaction rate for the scattering process 1 + 2 ^ JF into a set of final states J^ is related 
to the thermally averaged cross section times velocity, which can be calculated using the 
definition |l^ 



where 






dnr = g.^JiE,,T), (17) 



^We thank F. Feruglio on this point. 



Qi is the number of internal degrees of freedom for the particle species i and f{Ei, T) is the 
statistical distribution function for a particle with energy Ei in a thermal bath at temperature 

T. The factor VM01 = v (pi -^2)^ — ""^i^i/l-^i-^a) is known as the M0ller velocity (see 
Ref. [|12|) and a is the sum over all possible scattering channels of particles 1 and 2. 

In general equation (|TB|) can only be solved numerically, but if the particle species i 
is massless then it is possible to obtain analytic formulae for (|1^). Since the decoupling 
temperature of a particle species i will be calculated when T :» m, it will be a good 
approximation to use the analytic formulae for massless particles. For a cross section with 
a power law like behaviour and parametrised as a = J2i^is"'", where s is the Mandelstam 
variable one finds that for Re(nj) > —3 
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where Ci = 2^"^ — 1, ((2^+"^ — 1)^/18) for Bose-Einstein (Fermi-Dirac) statistics, ( is the 
Riemann zeta function and F is the Euler gamma function. If the cross section has a 
logarithmic dependence on s then for a = ao log s + ai with a boson and fermion in the 
initial state we obtain 



{crvM0i) = (To log T^ + 2.63327(To + ^i. 



(19) 



Let us first calculate the annihilation rate for the goldstino. In the limit that 
V^ ^ ^^3/2, '"^5, '"^p the dominant goldstino annihilation cross section is that into pho- 
tons (|l]). Using the fact that the equilibrium number density of relativisitic goldstinos is 



"g 



3C(3)/(27r^)T^ the goldstino annihilation rate is 



Tg = nG{avM0i) 






(20) 



Similarly, we can determine the interaction rate for the scalars 5* and P. Comparing the 
cross section ([T0|) (or (|13|)) with (|1^) we find that the scattering of S, P off the background 
thermal plasma (|T^) will be the dominant process determining the interaction rate. For the 
scattering off relativistic fermions with equilibrium number density Uf = 3C(3)/7r^T^ the 
scalar interaction rate becomes 
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and similarly for P. As discussed previously, nonzero cubic scalar couplings can affect the 
interaction rate at lower temperatures via the s-channel annihilation of scalars into two 
photons (p!B|). In this scenario with an equilibrium scalar number density ns = C{3)/n'^T^ 
the interaction rate is given by 
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(22) 



We will use this rate to make sure that the scalars do not come back into thermal equilibrium 
during nucleosynthesis. 
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FIG. 1. The Hubble expansion rate (solid line), Ygj (long-dashed line), Tss s-channel (S ex- 
change) (dashed line), Tss t-channel (^7 exchange) (dotted line) and Tq^ (dot-dashed line) as a 
function of temperature in the early universe. Notice that the goldstino decouples before the spin-0 
superpartners. 



Notice that for temperatures T ~ (9(100) MeV the scalar S and P interaction rate 
(pT|) is larger than the goldstino rate pO|). Consequently S and P will remain in thermal 
contact with the plasma much longer than the goldstinos. This can be clearly seen when 
the interaction rates are compared with the Hubble expansion rate of the universe H. If 
we express the energy density of the universe, p{T) in terms of the photon energy density 
P-y{T) = 7r^/15T^ then the Hubble expansion becomes 



H 



ISttG 



N 



-P 






(23) 



where Qp is the effective number of relativisitic degrees of freedom. In Figure |I] we have 
shown how the various interaction rates compare with the Hubble expansion rate for the 
representative values of gp = 65/4, M = 100 GeV, F = 3.8 x 10^ GeV^ and we have used 
M^ ~ ^/aT. It is clear from the figure that the goldstino will decouple at times much earlier 
than the spin-0 scalar superpartners and that the dominant scattering processes for S and 



P come from the scattering off background fermions (21). Also in the figure we can see the 
linear temperature dependence of (|2^) (with As = 0.01 GeV) and the possibility of F > if 
at lower temperatures. 

Let us now obtain the expressions for the decoupling temperature. In principle to de- 
termine the decoupling temperatures of S, P and G we would need to solve a set of coupled 
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Boltzinann equations. However, in practice since the interaction rates are significantly differ- 
ent between the scalar superpartners and the goldstino one finds that the coupled differential 
equations become separate equations controlled by ratios of T /H. In particular the Boltz- 
mann equations for S{P) separate when the dominant interaction rate is TsfiXpf) since 
the equations do not involve either of the other species. Thus it will be sufficient for our 
purposes to consider the equation F ~ i7 for the dominant interaction rate F. When the 
cubic scalar couplings are included we will also compare the interaction rate with H even 
though the equations may become coupled. Again this will be sufficient in order to obtain 
an order of magnitude estimate. 

When the cross section has power law behaviour the decoupling temperature may be 
easily obtained analytically. One finds for the goldstino decoupling temperature the result 

/ rp4 \ 1/5 

Tg ~ 0.77(7^^° -J- . (24) 

For the spin-0 superpartners the dominant scattering mode is V sj- In this case, an analytic 
expression for the decoupling temperature of the spin-0 superpartners of the goldstino can 
be obtained since for a plasmon mass My ~ a/oT the temperature dependence in the 
logarithmn cancels in T sf- Thus by requiring T sf ~ H the decoupling temperature is given 
by 
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Since in limit ^/s ^ ^3/2, ms, mp the scalars 5, P decouple later than the goldstino (see also 
the figure) we will use ( pSj) to set a lower bound on the scale of supersymmetry breaking F. 
Note that, without loss of generality we are also implicitly assuming that S and P decouple 
at the same time. 

When the cubic scalar couplings are included the cross section ([T5|) is again a power law 
and the analytical expression for the decoupling temperature is 

1/5 ,5 fM^Mp\ , , 

T ^ ^mg-^I^A\ (^-^ j • (26) 

We will need to make this temperature as low as possible in order that the fields S (and P) 
do not come back into thermal equilibrium during nucleosynthesis. 

It is also instructive to consider the limit msjfnp ^ ^/s ^ ^713/2 in which only the 
goldstino is assumed to be in thermal contact with the heat bath at temperatures T ~ 
0{MeV). As noted earlier, in this limit the goldstino annihilation cross section into photons 
changes to the form (^. Now the annihilation cross section into photons and fermions are 
comparable and qualitatively similar a (x s^. Including the electron, muon and 3 families of 
massless neutrinos for the fermions gives for the decoupling temperature 

T,^0.,3,r(^y\ (27) 



This equation will be used to show how much weaker the bound on the supersymmetry 
breaking scale F, becomes in the massive ms,mp case. 

In order to bound the decoupling temperature of a particle species i, we need to consider 
its effect on nucleosynthesis. This can be done by noting that the energy density of new 
massless particles is equivalent to an effective number AA^^ of additional doublet neutrinos 
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where fs = 8/7 for bosons, fF = ^ for fermions and Qi is the number of internal degrees 
of freedom of the particle species i |jl3|. If the scalars decouple at Ts^p = T,^ ~ (9(1 MeV) 



then according to ( P^ S and P will effectively behave as 1.14 extra neutrino doublets. If, 
however the scalars were to decouple during the epoch T^ < Ts^p < T^, where T^ is the 
muon decoupling temperature then the effective number of additional doublet neutrinos 
during nucleosynthesis due to S and P is 0.91. While neither of these bounds can be strictly 
ruled out [Q we will nevertheless suppose that S and P decouple before the muon decouples^ 
to quote a bound on the supersymmetry breaking scale F. 

As alluded to earlier there is no loss of generality if we calculate the bounds using the 
expressions for the scalar field S. The bounds are identical if we were to use the expressions 
for P, since to leading order the dominant cross sections have no dependence on ms or mp. 
Thus using the expression ( p5|) and requiring that Ts > T^ — (9(100 MeV) leads to the 
bound 

This strengthens the bound considered in M, where the effects of massless S and P during 
nucleosynthesis were not considered. For M ~ 100 GeV the bound on the supersymmetry 
breaking scale becomes F ^ 4 x 10^ GeV^. Using the relation F = v^^^3/2^p this leads 
to a bound on the gravitino mass 7713/2 ^ 1 eV. This is similar to the bound obtained by 
0] where the logarithmic term was neglected. One may have expected a large enhancement 
from the logarithmn factor in the massless limit due to the singularity, but since the plasmon 
mass is large, the enhancement in the cross section is only log a ~ —5 and thus impotent. 
Also if the mass of the scalars nis^p ^ 1 MeV then the bounds calculated here would be 
void because the cross section limits are no longer valid and a more detailed analysis in this 
intermediate mass range would need to be done. 

When the cubic scalar couplings, represented by As are included we need to make sure 
that they do not return to thermal equilibrium during nucleosynthesis. Requiring T ^ 0{eV) 
leads to a bound of 7723/2 ^ lO"^^^, which vanishes in the limit As = 0. We emphasise 
that while this bound is approximate, a more detailed analysis than that given here could 
provide interesting constraints on theories with nonzero cubic scalar couplings. 



^ This is somewhat plausible given that there may be other contributions to AN^, arising from 
massive neutrinos. 
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Finally, the gravitino mass bound coming from (^) in the limit ms,mp ^ ^/s ^ ^7^3/2 
is m3/2 ^ 4 X 10~^ eV and is independent of the gaugino mass M. As expected the bound 
becomes much weaker in this case and there are already stronger bounds from collider 
phenomenology in this limit [^]. 

IV. CONCLUSION 

We have shown that in the massless limit the scalar superpartners S, P of the goldstino 
couple much more strongly to the background thermal plasma than the goldstino. This 
means that during nucleosynthesis the scalars S and P may contribute as extra effective 
neutrino doublets, while the goldstinos have long since decoupled. Requiring that S and P 
decouple sufficiently early, so as not to upset current nucleosynthesis bounds, sets a lower 
bound on the scale of supersymmetry breaking which depends on the gaugino mass M. The 
bound on the supersymmetry breaking scale assuming M = 100 GeV is vF ^ 60 TeV 
which leads to a lower bound on the gravitino mass of 7713/2 ~ 1 eV. These bounds are 
much stronger than those considered in [Q, but complement the cosmological bounds in 
P]. Finally in the massive S,P case the nucleosynthesis bound is considerably weaker and 
similar to other recent cosmological bounds |p!5|,p!6|. 
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